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Abstract. We present a study of the one-particle spectral properties for a variety of models of Luttinger 
liquids with open boundaries. We first consider the Tomonaga-Luttinger model using bosonization. For 
weak interactions the boundary exponent of the power-law suppression of the spectral weight close to 
the chemical potential is dominated by a term linear in the interaction. This motivates us to study the 
spectral properties also within the Hartree-Fock approximation. It already gives power-law behavior and 
qualitative agreement with the exact spectral function. For the lattice model of spinless fermions and the 
Hubbard model we present numerically exact results obtained using the density-matrix renormalization- 
group algorithm. We show that many aspects of the behavior of the spectral function close to the boundary 
can again be understood within the Hartree-Fock approximation. For the repulsive Hubbard model with 
interaction U the spectral weight is enhanced in a large energy range around the chemical potential. At 
smaller energies a power-law suppression, as predicted by bosonization, sets in. We present an analytical 
discussion of the crossover and show that for small U it occurs at energies exponentially (in —1/U) close 
to the chemical potential, i.e. that bosonization only holds on exponentially small energy scales. We show 
that such a crossover can also be found in other models. 



PACS. 71.10.-w Theories and models of many electron systems - 71.10.Pm Fermions in reduced dimensions 



1 Introduction 

Theoretically it is well established that interacting fermi- 
ons in one spatial dimension do not obey Fermi liquid 
theoryjlj . The generic low-energy physics of one-dimensional 
(ID) metallic fermions with repulsive interaction can be 
described by Luttinger liquid (LL) theoryJl|,|,|,|,|,§. For 
various correlation functions it predicts asymptotic power- 
law behavior with exponents which, for spin rotational in- 
variant models, can be expressed in terms of a single pa- 
rameter K p . The Luttinger liquid parameter K p depends 
on details of the model considered, e.g. the interaction, fill- 
ing factor, and one-particle dispersion^, It has been de- 
termined for many different models of ID correlated elec- 
trons using either analytical or numerical techniques 
[jj . While the basic understanding of LL behavior emerged 
in the study of ID systems with periodic boundary con- 
ditions (PBC), the theoretical expectation that LL's with 
PBC including impurities scale to chains with open ends|)], 
|lC| led to several studies of models with hard walls, usually 
called "op en" (or "fixed") boundary conditions (OBC)|fll|, 
Il2|,[l3, 14, ill. One way to experimentally verify the pre- 



power-law suppression of the bulk spectral weight with an 
exponent 



dieted LL behavior is to probe the one-particle proper- 
ties using high resolution photoemission spectroscopy. For 
the spectral function p(u) entering the description of an- 
gular integrated photoemission and energies asymptoti- 
cally close to the chemical potential /i, LL theory predicts 



a=(K p + K; 1 -2)/(2z), 



(1.1) 



where z = 1 for spinless fermions and z — 2 for spin 1/2- 
fcrmions. 

For the translational invariant system the scattering 
processes which dominate the low-energy physics can be 
classified as forward, backward, and umklapp scattering ||. 
If the model parameters are such that backward or umk- 
lapp scattering become relevant in the renormalization- 
group (RG) sense both processes can drive the system into 
a gapped non-LL phase. In lattice models and for commen- 
surate filling factors umklapp scattering becomes relevant 
at a critical value of K p which dependents on the filling. 
For the lattice models considered here we will always chose 
the interaction and filling such that umklapp scattering re- 
mains irrelevant. Standard RG arguments can be usedj^] 
to show that for repulsive interactions the 2/c^-scattering 
part of the two-particle interaction (usually called u g\- 
interaction" ) scales to zero, where kp denotes the Fermi 
wave vector. Then the Tomonaga-Luttinger (TL) modelj^, 
H describes the generic low-energy LL physics. It only 
contains the scattering processes with small momentum 
transfer, which can be written as a quadratic form in 
bosonic density operators of left and right moving fermions. 
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This bosonization of the Hamiltonian is one of the ways to 
exactly solve the TL model with PBCfUjJ). For OBC we 
do not have momentum conservation and a general two- 
body interaction leads to a variety of different scattering 
vertices (see Sec. [2]). They depend on different combina- 
tions (differences and sums) of the four external quan- 
tum numbers and cannot simply be parameterized by the 
"momentum transfer" as for PBC. Only a few of the ver- 
tices can be written as quadratic forms in boson operators 
similar to the case of PBC@|. fn some of the previous 
publications on the open boundary problem pT|, p"2| it was 
tacitly assumed that RG arguments similar to the bulk 
case can be applied to show that the remaining vertices of 
the system without translational invariance scale to zero. 
Therefore a quadratic form in boson operators was used 
to describe the electron-electron interaction. Then it is 
straightforward to calculate correlation functions ]lj|,|l2| ■ 
In Sec. ^ we present a detailed discussion of bosonization 
for the case of OBC taking all the scattering processes 
into account. We explicitly demonstrate that the above 
assumption is justified for an interaction which is long 
range in real space, considered by Tomonaga for PBCQ. 
In Refs. |Tl| and [[12| it was shown that the local spec- 
tral density p(x,co) near the end points of a ID chain is 
modified compared to the bulk density. The algebraic be- 
havior of the spectral density with frequency u close to the 
chemical potential was found to be governed by a bound- 
ary exponent 



CtB 



(K- 1 - l)/z, 



(1.2) 



which, for repulsive interaction (K p < 1), is larger than 
the bulk exponent a. 

Systems of ID correlated electrons can be viewed as 
being at a quantum critical point H. The occurrence of 
nonuniversal (critical) exponents can then be traced to 
the fact that the low-energy physics is governed by a line 
of fixed points. From the theory of critical phenomena it is 
known that critical exponents at a boundary usually differ 
from their bulk counterparts. Having this in mind the ob- 
served difference between olb a nd a is not surprising 

To investigate whether Eq. (1.2) also characterizes the 
spectral function for models with OBC and a short range 
interaction we have calculated the spectral weight at the 
boundary site and the chemical potential using the density- 
matrix renormalization-group (DMRG) method0 for the 
lattice model of spinless fermions with nearest neighbor 
interaction U and the ID Hubbard model with onsite in- 
teraction U . In Sees. || and |^ we show that the numerically 
exact data are indeed consistent with Eq. ([L2]), although 
for the Hubbard model only at energies extremely close to 
p. 

For small interactions and the models considered K p — 
1 can be expanded to give a leading behavior which is 
linear in the interaction Jl|,@,||]. For a this gives a lead- 
ing term which is quadratic in the interaction. In second 
order perturbation theory for the self-energy the nonana- 
lytic power-law behavior appears as a logarithmic diver- 
gence In \uj\ with a prefactor which is of second order in 
the interaction^]. In contrast olb has a contribution lin- 



ear in the interaction. Thus signs of the nonanalytic be- 
havior of p(x,co) can already be obtained using the (non- 
self-consistent) Hartree-Fock (HF) approximation for the 
self-energy. The least we can expect to find within the HF 
approximation is a logarithmic divergence of the above 
form with a prefactor linear in the interaction. In Sees. ^, 
[|, and |^ we study the spectral function for the three mod- 
els (TL model, spinless fermions, Hubbard model) consid- 
ered within the HF approximation. By comparison with 
the exact results (bosonization and DMRG) we show that 
many aspects of the low-energy behavior of p(x, uj) can 
be understood within the HF approximation. Surprisingly 
p HF (x, ui) already gives power-law behavior. This is a very 
interesting observation as the HF approximation for the 
case of PBC does not capture any of the LL features. 

Neither from general LL theory nor from the theo- 
retical and numerical analysis of microscopic models, e.g. 
the ID Hubbard model, much is known about the energy 
range A over which the power-law behavior in the spectral 
function can be observed. Obviously a knowledge of A is 
essential for a meaningful comparison of theoretical and 
experimental spectra. The analytical techniques used, e.g. 
bosonization and the Bethc Ansatz in combination with 
boundary conformal field theoryjl]], only provide the ex- 
ponent which characterizes the spectral weight at energies 
asymptotically close to the chemical potential. This holds 
for both PBC and OBC. Quantum Monte-Carlo calcula- 
tions of the spectral weight for the ID Hubbard model 
with PBC indicate that A is small[|l8|], which implies that 
chains of many lattice sites are required to observe the 
suppression in numerical calculations. In Sec. ^ we show 
that for the repulsive Hubbard model with OBC the spec- 
tral weight close to the boundary is enhanced in a large 
energy range around the chemical potential. The power- 
law suppression, predicted by bosonization, only occurs 
after a crossover at energies A which for small U are 
exponentially (in —1/U) close to the chemical potential. 
For the particular case considered this demonstrates that 
the bosonization (and conformal field theory) result only 
holds on a very small energy scale. We present an analyti- 
cal discussion of the crossover behavior within an effective 
model. From perturbation theory we expect that this kind 
of crossover occurs in all models with a bare total back- 
ward scattering amplitude, i.e. z times the 2fcp component 
of the interaction, which is larger than the forward scat- 
tering. 

Several attempts have been made to experimentally 
demonstrate LL behavior in a variety of systems which 
behave as quasi one-dimensional conductors using angular 
integrated and angular resolved photoemission spectrosco- 
py |lj],^0],^l],|2^,^3|. Unfortunately all the measurements 
are plagued by various subtle problems [Ts[| and their in- 
terpretation has been questioned [ ^9||22]| . Until three years 
ago the experimental spectra were compared to the pre- 
dictions of bulk LL theory. Using the results obtained 
from bosonization and conformal field theory very recently 
many authors argued that the theoretical picture of a 
chain which, by impurities, is cut into several disconnected 
pieces gives a better agreement between the spectra ob- 
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served in photoemission experiments and theoretical spectra| 
[l^, ^9| , |20[|2^| . In Sec. fj] we will discuss this issue based on 
our results for the spectral function of the Hubbard model 
and the expectation that a crossover behavior similar to 
the one observed for the Hubbard model occurs also in 
other microscopic lattice models. Our results demonstrate 
that features which occur at energies not captured by 
bosonization and conformal field theory might dominate 
spectra broadened by finite temperatures and experimen- 
tal resolution. Thus a reliable interpretation of photoemis- 
sion data requires further theoretical investigations. 



2 Open boundaries and bosonization 

In this section we present a detailed discussion of bosoniza- 
tion for a continuum model of length L with OBC and a 
two-body interaction with a spatial range R = l/q c - We 
discuss the subtleties which in previous approaches have 
only partly been considered |ll[|l2] ]. 

The one-particle eigenstates of the noninteracting sys- 
tem are given by 
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tp n (x) = \/2/L sin (k„x), 



(2.1) 



with k n = mr/L, n 6 N. Note that in contrast to the 
translational invariant case the k n do not have the mean- 
ing of momenta. A ID system with PBC has two Fermi 
points ±fcj?. In contrast here we only have one Fermi point 
given by kp = npir/L, where znpp4[ denotes the number 
of electrons in the system (z = 1 for spinless fermions and 
z = 2 for spin 1/2- fermions). The interaction between two 
particles of spin species s and s' is characterized by a two- 
body potential V s>s i{x — x') which leads to a contribution 
to the Hamiltonian given by 

V= ± J2 [ L dx f L dx'i,t(x)^ sl (x') 

1 S,8' J ° J ° 

xV StS '(x - x')ip a ,(x')ip s {x) 
= iV / dx I dx'p s (x)V s>s '(x-x')p s >{x') 

1 S,8' J ° J ° 

s 

with the field operators ,(x) — J2ri=i ( Pn(x)an} s , the 
density operators p s (x) = ipl(x)ip s (x), and the particle 

number operators N s = J dxp s (x) . We express the inter- 
action in terms of the creation and annihilation operators 

tJn,l of the eigenstates tp n with spin s. After rearranging 
the terms in a way which simplifies the bosonization dis- 
cussed later it reads 

^ = 2 V nmrim' a li,s a n' ,s a m,s' a m' ,s' 



n,m 



nmnm 



with the matrix elements 

L 



dx 



o 



5 S , S ,V S , S (0)N S 



dx'ip* n {x)ip n ,{x) 



(2.3) 



xV s , s >(x - x')tp* m (x')ip m ,{x') 



(2.4) 

and the occupation number operators n„ jS = a], s a n s . If 
we express products of the sine functions ip n in terms of 
cosine functions they read 



nmn m- 

-F, 



[-^s.s r {kn 

{k n k n > , k m 



kn' ? km 



Fs.s' {kn 



kn' ; k m 



km' ) 

km') ~~ -^s.s' {kn H~ k n f ikjyi &m') 

■k m ,)]/L, (2.5) 



with 



F StS i(q,q') = — I dx dx' cos (qx)V s , s *{x - x') 
L Jo Jo 



x cos (q'x') 
1 

4L 

xV SlS '(x - x 1 ). 



dx dx 



(2.6) 



The area D over which we have to integrate in Eq. (2.6) is 
given by the two hatched squares shown in Fig. 2.1 each of 
size L 2 . Using the fact that V s , a >(x) was assumed to have 
a range R, F s s > (q, q') can partly be expressed in terms of 

the Fourier transform V s ,s' (q) = J^ L dxV StS ' (x) exp (— iqx) 
if we instead inte grate over the rectangle of width 2R indi- 
cated in Fig. 2.1. This can be achieved by adding and sub- 
sequently subtracting the integral over the six black tri- 
angles. If we denote their contribution by —g SfS i(q,q')/L, 
we obtain 

F s ,s>(q,q') = V s ,s>(q)/2 (<W + Sq-q') + g s ,s>(q,q')/L. 

(2-7) 

In the case of a spin independent exponential interaction 
V(x) = Vb/(2i£) exp (— \x\/R), g(q,q') is proportional to 
V(q)V(q')/V . ' ' 

For comparison we give a brief discussion of the in- 
teraction part of the Hamiltonian for the ID electron gas 
with PBC. The noninteracting eigenstates are plane waves 
4>n{x) = exp (ik n x)/y/~L with k n = 2-KnjL and n € Z. In 

terms of the creation and annihilation operators &n,s of 
these st ates a two-particle interaction gives a term similar 
to Eq. (2.3) with the matrix elements 



~s,s 
nrnn' m 



Fs : s' {krt 



h , k 



km' ) 



and 



F StS ,(q,q') = V SlS ,(q) <W/2. 



(2.9) 



For PBC the matrix elements thus only depend on the 
momentum transfer k„ — Av = km> — k m between the 
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Fig. 2.1. Area of integration for the integral in Eq. (2.6). For 
details see the text. 



two scattering particles. The part of the interaction with 
small momentum transfer q <C fc^ can be written as a 
bilinear form in the spin and charge density operators of 
left and right moving fermions @, [||lj . There are two ways 
of treating the backscattering processes with momentum 
transfer 2fcp. One can show that for repulsive interactions 
backscattering is irrelevant in the RG sense Q. If one is 
only interested in the "critical" low-energy properties of 
the model one can thus neglect backscattering from the 
beginning. Then the physical two-body potential has to be 
replaced by effective coupling constants, as it is changed 
by the initial flow of the backscattering and one ends up 
with an effective low-energy model Q. The other way to 
proceed is to assume that the interaction V s>s i (q) is cut 
off at a momentum q c = l/R <C kp, i.e. is long range in 
real space and only consider a subspace of the Fock space 
Tt with no holes deep in the Fermi see and no particles 
in one-particle states with energies much higher then the 
Fermi energy. This is the idea originally considered by 
TomonagaB. On this subspace the properly normalized 
density operators of left and right moving fermions obey 
bosonic commutation relations |lq| . For a one-particle dis- 
persion which is linearized around the two Fermi points 
also the kinetic energy is quadratic in the densities fig] and 
the entire Hamiltonian can be written as a Hamiltonian 
of noninteracti ng b osons. 



From Eqs. ( |2.5| ) and (2.7) it is obvious that the situ- 
ation is more complex for the case of OBC. Expressed in 
the basis of the noninteracting eigenstates the interaction 
contains a variety of scattering processes. Only a few of 
them can be written bilinearly in the density operators 
(m 6Z,m/0) 



Pm,a = (/?m,T _ Pm,l) / V2 



(2.12) 



can be used. These operators are defined in analogy to 
the case of PBC. A RG study which includes all the scat- 
tering vertices is still missing and we will thus follow the 
idea of Tomonaga to investigate whether this leads to a 
Hamiltonian bilinear in the densities similar to PBC. Thus 
we assume that the interaction V SiS > (x) is long range in 
real space and only consider the low-energy subspace Tt- 
For simplicity we will furthermore focus on the case of a 
spin independent interaction V s>s i(x) — V(x). Acting on 
states IV't) € Tt the interaction term V simplifies because 
a n s a n',s \iPt) = for < \n — n'\ < n c , unless k n i is close 
to kp. As, with our assumption about the range of the 
interaction, F(q, q') is only nonzero if both arguments are 
smaller than q c we can rep lace v nmn i m i in the first term on 
the right hand side of Eq. ( |2.3| ) by F(k n — k n / , k m — k m >)/L. 
With similar arguments for the other contributions the in- 
teraction term on Tt reads 

1 v^' 



yJ2 F P^> ^p(Pn,p + P-n,, 



n>0 

fx 



n>0 m>0 



^F p (0, 0)iV p 2 - ^ E F P& 2k ^P 

m>0 

^-V p (0)N p + j £ F„(2A: m ,2fc m /),(2.13) 



with 



N, 

F P (q,q 
v P (Q) 



m.m'X) 



p -(Ni+N L )/V2, 
" = zF s>s ,(q,q'), 



zV s , s ,(0). 



(2.14) 
(2.15) 
(2.16) 



The prime at the sum in Eq. (2.13) indicates that the 
terms with n = or n! = are excluded. For z — 1 the 
above formulas give the corresponding expressions for the 



case of spinless fermions, if one puts p 



and N n 



N. Co mpared to the analogous expression for PBC|16| Eq. 
(2.13) contains three modifications: Generically F p (q n , q n i) 
is nonvanishing for all q n and q n i smaller than q c which 
leads to a co upling between all the p„ jP in the first line 
of Eq. ()2 13|)|j rj|. For PBC only p n , p with n = -ri are 



n— max{ 1,1— m} 



n.s n+m,s ' 



(2.10) 



which obey p^ s - 
density operators 



coupled. OBC furthermore lead to a term which couples 
the particle number and the density operators and a term 
linear in the p w ,p ]25|] • As we will discuss in Sec. ||the lat- 
ter leads to a nontrivial contribution to the Hartree term 
in the self-energy. For the discussion of the one-particle 

Alternatively charge and spin properties the constant as well as terms linear in N p can 

be neglected as the latter only lead to a renormalization 
of the chemical potential|jl^|. If we linearize the disper- 

= Pn.v/\/n for 



(Pm,T + Pm,i)/^ 



(2.11) sion around kp and for n > define b Tl 
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v — p, er the kinetic energy on Tt can, up to particle num- 
ber contributions, be replaced by 

H =V F jJ2 n ( b lp b n, P + «4,<A,<r) (2-17) 

with the Fermi velocity vp. On Tt the b n p ^ ij obey bosonic 
commutation relations and charge and spin degrees of free- 
dom do commute. Note that in H = Hq + V similar to 
PBC spin and charge degrees of freedom are decoupled. 
Formally a Hamiltonian H of coupled and shifted har- 
monic oscillators can be diagonalized using a Bogoljubov 
transformation. For a general interaction the transforma- 
tion cannot be given analytically, but e.g. for the above 
mentioned interaction V(x) = Vq/(2R) exp (— ex- 
ploiting the fact that the corresponding g(q, q') is separa- 
ble, it can analytically be shown that neither the coupling 
of bosons with \q n \ ^ \q n '\ given by g(q n , <2Vi') nor the terms 
linear in the bosons do change the "critical" one-particle 
properties of the model, i.e. the ex pon ents of the asymp- 
totic decay of correlation functions |25| . As long as we are 
only interested in the "critical" behavior we can thus work 
with the simplified Hamiltonian 



H = [v F ^bl p b n p + ^V(q n ) (bl p + b n p ) 

•j^J2 nb ^ b n 



n>0 
-Vp- 



(2.18) 



n>0 



In the following we will denote the model related to this 
Hamiltonian the TL mo del, ev en though we have dropped 
terms compared to Eq. ( 2.13 ) as desc ribed above. To ob- 
tain the spinless model from Eq. ( p. 18 ) one has to drop the 
second line and set z = 1. The TL model with OBC con- 
sists of an independent system of selfcoupled oscillators. 
This is in contrast to the PBC case, where the modes n 
and —n are cou pled by the interaction term. Apart from 
a constant Eq. ( p. 18 ) can be brought into the formpd|,^6| 



H 



n>0 



v F^nbl ia b niJ 



(2.19) 



with uj n = k n yl + zV (k n ) / (jrv f) and bosonic operators 
a n , P given by a linear combination of b' n „ and b n „. 

The one-particle spectral function p < (x, uj) relevant for 
photoemission can be determined from the Green's func- 
tion G<(x,x,t) = -i(ipl(x,0)il) s (x,t)) using|2|l 



p < (x,u) 



dt 
2^ 



e iuit iG < (x,x,t). 



(2.20) 



Here (. . .) denotes the ground state expectation value. For 
the TL model with OBC G < (x, x, t) can be calculated 
using the bosonization of the fermionic field operator 



L ^ 

=i 



sin (k n x) a n 



<2L 



E 

n=l 



(2.21) 



Eq. ( 2.21 ) cannot directly be expressed in the bosons bn, P 

and bn,a- To achieve this we have to add one-particle states 
with quantum numbers k n , n < 1 to the Hilbert space. 
They are assumed to be filled in the ground state and 
do not modify the low-energy physics of the model. The 
auxiliary field operator 



oo 



(2.22) 



can then be written in the boson operators similar to 
PBCjl6) and G < (x,x,t) is given by 

G < (x,x,t) - [^1^,0)^^,^) + ^t(~x,0)^ s (-x,t)) 

^J(x,0)^(-x,*)\ - Ut(-x,0)$.(x,tj)] 12. (2.23) 



The expectation values in Eq. ( |2.23| ) can be calculated 
using bosonization of the auxiliary field operator ]l6|. For 
a fixed position x the leading behavior of G < (x,x,t) at 
large times isO] 



G<{x,x,t) 



with 



a B = (k; 1 - l)/z 
and the LL parameter 

-1/2 



1 



zV(0) 
ttvf 



(2.24) 
(2.25) 

(2.26) 



of the TL model. The nonanalytic behavior of the spectral 
function close to the chemical potential is thus given by 



p<(x,u)~\u\ aB G(-w). 



(2.27) 



Without explicitly demonstrating that the TL model 
is the effective low-energy model (fixed point model ) for 
all mo dels of LL's with OBC the results of Eq s. (|2.25 ) and 
( [2.27D have been assumed to hold for all LL'sjOj]!^ This 
implies that the boundary exponent as can be expressed 
in terms of the bulk LL parameter K n . The generaliza- 
tion found some confirmation in Ref. [13| where methods 
of boundary conformal field theory were used to calculate 
as for Bethe ansatz solvable models. In Sec. |^ we will 
present numerically exact results for the spectral weight of 
the lattice model of spinless fermions with nea rest n eigh- 
bor interaction, which are consistent with Eq. (2.27) over 
a fairly large energy range. Our results for the Hubbard 
model presented in Sec. H demonstrate that in this model 
the spectral function displays a richer structure. Only for 
energies A exponentially (in —1/U) close to fi, i.e. expo- 
nentially large syste m sizes, the small U data are consis- 
tent with Eq. Q2.271 ). 

Within the TL model and for fixed interaction strength 
the x dependence of the energy range A over which the 
asymptotic power-law suppression with the boundary ex- 
ponent as can be observed is given by 
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A k, vp min{l/x, 1/-R}- For x > R and vp/x < uj < vp / R 
the bulk "critical" behavi or w ith a power-law suppression 
given by a defined in Eq. ( |l.l| ) is recovered. For u> > vp/R 
"nonuniversal" features dominate the spectrum. Unfor- 
tunately the energy scales obtained from the TL model 
cannot be trusted if one is interested in more realistic 
continuum or lattice models with a nonlinear one-particle 
dispersion and scattering processes given by a general in- 
teraction (not necessarily long range in real space). The 
best bosonization can provide for such models is the ex- 
ponent of the nonanalytic behavior at energies asymptot- 
ically close to ji. 

For small V(Q) it follows from Eq. fl2.26|) that 



ith 



zV(0) 

2ltVF 



O 



no) 

TTVF 



and thus from Eq. ( |2.25| ) that 



a B 



no) 

2nvp 



O 



ttvf 



(2.28) 



(2.29) 



This has to be contras ted to the small V(Q) behavior of the 
bulk exponent a Eq. (1.1) which is quadratic in the inter- 
action. Thus signs of the nonanalytic behavior of p[x 1 ui) 
can already be obtained using the HF self-energy, which 
will be analysed in detail in the next section. 



3 Hartree-Fock self-energy 

In this section we discuss the non-self-consistent HF ap- 
proximation for the self-en ergy (S^ p ) n ^ n i , for the contin- 
uum model defined by Eq. (2/2) and a one-particle disper- 
sion e(k). In contrast to the end of the last section we do 
not restrict our discussion to the subspace Tt and thus 
con side r all the scattering processes given by the Eqs. ( |2.3| ) 
to (2.7). We assume that the interaction is spin indepen- 
dent. For PBC (S^ F ) n . n i oc 6 n>n > because of momentum 
conservation and the HF approximation only leads to fi- 
nite shifts in p, and Vp. It does not capture any of the pe- 
culi ar pro perties of LL's. As al read y mentioned following 
Eq. ( 2.17 ) g s , s '(g, q')/L in Eq. (2/7) does not contribute to 
the "critical" LL properties of the continuum model with 
OBC. Later we will confirm this using perturbation the- 
ory. We will thu s ne glect this term and only consider the 
5 ter ms in Eq. (2.7). With the matrix elements given by 
Eq. (2J3) this leads to 



n" — 1 s' 

n,n' I Sp, - — ^V(k n - k ni ) + V(k n + k ni ) \ 

I 711=1 J 



n + n' 



f(m) = 



1 for m G {1, 2, . . . , uf} 
otherwise 



(3.2) 



and 5/i = zV(0)np/ L. Due to the parity symmetry with 
respect to the middle of the box only matrix elements with 
even n ± n' are nonvanishing. Note that the self-energy is 
to independent and real but has a nontrivial matrix struc- 
ture in the quantum numbers n and n' due to the broken 
transnational invariance. From the self-energy the retarded 
Green's function follows by a matrix inversion pcj 

[G(f>)Un> = " ^) + iO} 1 - S}-%, , (3.3) 

with the unity matrix 1 and £(fc) = e(k) — /i. The local 
spectral function is then given by 

^ oo 

IT ' 



p < (x,uj) defined in Eq. ( | 2.20 [) is related to the total spec- 
tral function of Eq. (|3.4|) by p < (x,uj) — p(x,u)0(uj). To 
specify the model we have to chose a dispersion e(k) which 
is assumed to be sufficiently smooth. In the thermody- 
namic limit the noninteracting spectral weight is given 

byH 



p°(x,u) 



l 



ds(k) 



dk 



[1 - cos [2kx)\ 



(3.5) 



where ^ (w) denotes the function inverting £(fc). We did 

not suc ceed in analytically inverting the matrix [G HF (u))] 1 
in Eq. (3.3) which, for a general interaction, has a nontriv- 
ial structure. To gain a first insight in the behavior of the 
spectral weight close to /i the matrix can be inverted per- 
turbatively to lowest order in V 



p HF (x,c)=p°(x,c) 1 + 



V(0) V(2k F ) 
2nvp 



2lTVf 



x In 



VFk} 



O V 2 



(3.6) 



Note that the logarithmic divergence in p HF (x,uj) is not 
due to a singular frequency behavior of Z" HF , but emerges 
in the perturbative approach to the matrix inversion. This 
has to be contrasted to PBC where the first indication 
of a break down of perturbation theory can be found in 
second order. The second order self-energy S 1 - 2 ^ displays a 
logarithmic divergence V 2 In \uj\ leading to the same kind 
of divergence in the spectral function|5|,|l6||. The diagram 
responsible for the second order divergence for PBC does 
also lead to a logarithmic divergence 



,V 2 lnU 



(3.7) 



for OBC. Thus logarithmic terms found in perturbation 
theory for p(x, ui) in the case of OBC can have two differ- 
ent origins. 
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The leading logarithmic behavior in /? HF (x, to) Eq. (3.6) 
comes from the term in Z' HF which is proportional to 
f([n + n']/2). It gives a sharp step in the self-energy ma- 
trix [see Eq. ( |3.2|) 1 which crosses the diagonal at (nj?,nf). 
The height of the step determines the prefactor of the log- 
arithmic divergence. Contributions of g{q,q') neglected in 
Eq. (3.1) are continous and thus do not modify the non- 
analytic behavior, consistent with the observation made 
in the last section. In Sec. [(] we will calculate the HF self- 
energy for the Hubbard model with OBC and present an 
analytical discussion of the resulting spectral weight in an 
effective model. Sharp steps in the HF self-energy matrix 
will also play a prominent role in this discussion. 

To investigate how the leading logarithmic divergence 
of p tLF (x,uj) Eq. (3.6) is modified by higher orders in V, 

we numerically invert [G hf (cj)] for different models of 
finite size in Sees. |], || and It turns out that already 
the HF approximation leads to a power-law. We will fur- 
thermore compare the spectral weight in the HF approx- 
imation to exact results obtained from bosonization and 
DMRG. 



4 Tomonaga- Luttinger model 

The spectral funct ion f or the TL model of finite length 
L defined by Eq. ( 2.18 ) can be calculated exactly using 
bosonization and a recursive method introduced in Ref. 
p8| for the case of PBC. For the momentum integrated 
spectral function of the TL model the spin only leads to 
factors of two and we will thus c onsid er the spinless TL 
model. The last two terms in Eq. ( [2.23 ) give contributions 
to p(x, to) which for energies close to the chemical poten- 
tial are proportional to cos {2kpx). Such a term is alr ead y 
present in the noninteracting spectral function Eq. (3.5). 
It drops out if the spectral function is spatially averaged 
over a small length. If one is interested in the comparison 
to experimental spectra this assumption is justified as pho- 
toemission measurements will auto mat ically average over 
some small spatial range. In Fig. 4.1 only the nonoscil- 
latory part p no (iE,a;) of p(x,tu) is shown (solid line) as a 
function of uj. The low-energy behavior is independent of 
the shape of the interaction and we have chosen the sim- 
ple form V{q) = U0(q c — \q\). For a finite system p(x,u)) 
is given by a sum of S peaks. In Fig. 



1.1 



the energies with 
nonvanishing weight are equidistant with a level spacing 
oc 1/L and we have connected them to a continuous line. 
A power-law behavior of the weight close to p can only be 
found in the thermodynamic limit, but for large system 
sizes the weight of the 8 peaks resemble s th e power-law 
obtained for L — > oo. The data in Fig 4.1 display the 
suppression of spectral weight close to p which is much 
stronger then the one observed for PBC. For compari son 
PBC data for the same parameters are shown in Fig. 4.1 
as a dashed-dotted line. 

The HF s elf-en ergy f or th e TL model differs for the two 
choices Eqs. ( 2.13| ) and ( 2. 15 ) for the interaction. The term 
proportional to V(k n — k n ')f(\n — n'\/2) in Eq. (3.1) is due 
to the interaction term linear in the bosons neglected in 



1.0 



0.9 



> 



0.8 



-1.0 



-0.8 



exact, OBC 
Hartree-Fock, OBC 
exact, PBC 



-0.6 -0.4 
co/(q v) 



-0.2 0.0 



Fig. 4.1. Spectral weight near the boundary (x = 4i?) for the 
spinless TL model with U/(2ttvf) = 0.05263, n c = q c L/n = 
160, and uf = 4n c . 



Eq. (2.18). As it does not contribute to the singular part of 
the HF spectral funct ion w e treat the HF approximation 
to the TL model Eq. fl2.18p. This leads to 



l ni=l ) 



1 - / k„ — k„i \ I n + n' 
2L V M — 



(4.1) 



The first term in Eq. ( |4.l| ) is the HF shift of the chem- 
ical potential and the second term leads to a renormal- 
ized Fermi velocity i>hf — vf + V(0)/(2w). The resulting 
matrix [G HF ] which has to be inverted is sketched in 
Fig. [L^. If we consider the box shaped potential V(q) = 
UO(q c — \q\) the last term in Eq. (4.1) doe s contribute 
as a constant in the hatched area. In Fig. 



4.1 



HF 

rno 



(x, lu) is 

shown as the dashed line. For the parameters chosen (weak 
interaction) the HF result agrees quantitatively with the 
exact soluti on. F or the TL model V(2kp) = 0, which leads 
to [see Eq. (O)] 



2-kvf 



In 



to 



Vpk[ 



(4.2) 




Fig. 4.2. [G HF ] for the TL model. For details see the text. 
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For repulsive interactions (V(0) > 0) the prefactor of 
the logarithm is positive and the perturbative expression 
indicates a suppression of the weight. To further inves- 
tigate the behavior of p(x,oj) close to the chemical po- 
tential we have studied the spectral weight at the chem- 
ical potential and position x, denoted wo{x,np]V), for 
a given kp as a function of 1/np (oc 1/L). The ratio 
wq(x 1 np;V)/wo(x,np;0) displays the same kind of power- 
law behavior as does p(x,tu) as a function of to. Smaller 
systems are sufficient to find power-law behavior in 
wq{x, np] V)/wq{x, np] 0) compared to p{x,Lo). Fig. 4.3 
shows a log-log plot of wq(x, rip; V)/wq(x, rip; 0) for the 
exact solution and the HF approximation. A powe r-law 
fit of the s ymbols gives the expected exponent Eqs. ( 2.25 ) 
and ( 2.26 ) with high accuracy. Surprisingly also the HF 
approximation displays power-law behavior, which shows 
that the leading logarithmic divergence of p HF found in 
Eq. (4.2) can be resummed to give a power-law. A de- 
tailed study shows that a§ F = V (0) / (2ttvhf) and thus 
a% F and as do agree up to leading order in V(0)/(2nvp ). 
Quantitative agreement between exact results and HF can 
be reached for V(0)/(2irv F ) < 1. 



a o.8 



0.6 



• bosonization 
Hartree-Fock 

■ bosonization 
Hartree-Fock 



10" 



10" 



10~ 



1/n c 



Fig. 4.3. Spectral weight at /i and close to the boundary. The 
circles show the exact results for the spinless TL model for 
x = 3R, U/(2tivf) = 0.1 and kp — 4q c . The long dashed line 
presents the HF approximation. The squares show the exact 
results for x = 37? and U/(2hvf) = 0.05 and the dashed line 
the corresponding HF approximation. 



5 Lattice model of spinless fermions 

Next we consider the lattice model of spinless fermions 
with N lattice sites, lattice constant a = 1, hopping matrix 
element t = 1, nearest neighbor interaction U, and OBC 



H 



N-l 



: jS'+J 



N-l 

u n J n J+i- 



(5T) 



c^p denotes the creation (annihilation) operator at site j 
and rij = CjCj. For U = the eigenstates of H are given 



by Eq. ( |2T| ) with L ->■ 7V+1, x -> j and k n = mr/(N+l), 
n E {1,2,..., N}. Similar to PBC the one-particle disper- 
sion is e(k) = — 2cos(fc). In contrast to the TL model 
the interaction is of short range (in real space) an d it i s 
not obvious whether the bosonization result Eq. ( p.2li ) 



as = K^ 1 — 1 holds. As for PBC the interacting model 
with OBC can be solved exactly by the Bethe ansatz, but 
similar to PBC not much about correlation functions can 
be learned directly from the solution. Information about 
boundary exponents can be obtained if conformal invari- 
ance is assumed ]l3|,^5[. 

Here we discuss the local spectral function at site j = 
1. We have performed a DMRG study 10 for chains of up 
to N = 512 sites calculating matrix elements wo(riF; U) = 
K.B^"^!^^)! 2 , i. e. the spectral weight at the chem- 
ical potential and the boundary site (see Sec. ||). \Eq f ) 
denotes the exact n_F-particle ground state. Results for 
wo(np] U)/wo(riF] 0) as a function of 1/./V (instead of 
I/uf), two different filling factors n/ = np/N, and U 

K p (U,rif) for these 



are shown as the symbols in Fig. 5.1 



parameters can e.g. be found in Ref. 
0.5) = 0.9691 and K p (U = l,n f 



K p (U 
■■ 0.25) 



0.1, n/ = 
= 0.8447. 



For large N the numerical data nicely follow the solid 
lines, which are proportional to power-laws with expo- 
nent as(U = l,ri/ = 0.25) = 0.1838 respectively as{U = 
0.1,71/ = 0.5) = 0.0319. The numerical error of the DMRG 
data is smaller than the size of the symbols. We can thus 
conclude that the spectral weight close to the boundary 
and the chemical potential shows a suppression with a 
power-law and exponent as- This is consistent with the 
prediction of bosonization and results presented in Ref. 
[[f| obtained from Bethe ansatz and boundary conformal 
field theory. 

In contrast to the TL model the HF problem for spin- 
less fermions is numerically best studied in the site rep- 
resentation. The HF or mean-field Hamiltonian is given 

by 



N-l 



H 



HF 



E (t(j)c]c 3+1 +t*(j)c} +l( 
i=i 

N-l 

J2[U + u(j)} nj 



(5.2) 



with a renormalized hopping 

t(j) = l + U{c] +lCj ) 



= 1 



onsitc energies 



,7 + J 

2U 



n f 



N+ 1 ^ 



(k n j)sm(k n [j + l}), (5.3) 



(n 2 ) forj = l 

= U { (n,-_i) + (T»i+i)o for j e {2, 3, . . . , iV - 1} 
(n N -i) for j = N, 



and 



(%)o 



N 



— ^sin 2 (k n j), 



(5.4) 
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where 



denotes the noninteracting ground state ex- 



pectation value. Due to the OBC the site occupation 



depe nds on j s how ing Fricdel oscillations. The sums in 
Eqs. (5J3) and (5T) can be performed analytically and in 
order to obtain the site diagonal Green's function 
[G HF (a;)] and thus the spectral function p^ F (uj), one 
numerically has to determine the eigenvalues and eigen- 
vectors of a tridiagonal matrix. This can easil y be done 
for systems of up to N — 10 6 lattice sites. Fig. 5T shows 
HF data for wo(np; U)/wo(jif] 0) and the same param- 
eters as above. Similar to the TL model already the HF 
approximation displays power-law behavior with exponent 
a B F — Kfif/(27TUHp) which has the same form as for the 
TL model if one replaces V^(0) by the effective interaction 
V cS = V(0) - V(2k F ) = 2U[1 - cos(2fc F )]. This replace- 
ment is familiar from the mapping to the "g-ology" model 
for PBC(^]: For spinless fermions forward gi and backward 
gi scattering are indistinguishable which leads to an effec- 
tive coupling gi — g\ . Again a b and o^ F agree up to order 
V c s /(2irvp) and both curves show quantitative agreement 
for V eS /(2irv F ) < 1. 



0.8 



~ 0.6 



0.4 




• U 

- u 
■ u 

u 



0.1, n f =0.5, DMRG 

0. 1, n f =0.5, HF 

1, n,=0.25, DMRG 
1, n,=0.25, HF 

(1/N)° B 



1(T 



10" 



10" 



1/N 



Fig. 5.1. Spectral weight at /i and lattice site 1 for the lat- 
tice model of spinless fermions. The symbols show the DMRG 
data and the dashed lines the corresponding HF results. For 
comparison the solid lines are power-laws with exponent ats — 



K- 



1. The parameters are given in the legend. 



6 Hubbard model 

The main part of our earlier publication p9| on LL's with 
boundaries is devoted to a numerical investigation of the 
spectral properties of the ID repulsive Hubbard model 

N-l 

s j=l 
AT-1 

3=1 



Similar to the lattice model of spinless fermions Eq. (5.1) 
we set the lattice constant a and the hopping matrix ele- 
ment t equal to unity. Here we will give a short summary 
of the numerical results obtained in Ref. j2{| , discuss new 
data sets, and then present an analytical investigation of 
the crossover behavior of the spectral weight within an 
effective low-energy model. 



6.1 Numerical results: DMRG 

As for the lattice model of spinless fermions we have cal- 
culated the spectral weight at the boundary site and the 
chemical potential as a function of the number of lattice 
sites N using the DMRG algorithm. We were able to ob- 
tain results for up to N — 256 sites. Figs. 6.1 and 3.2 
show data for wo(nF\U)/wo(riF',0), quarter filling nj = 
2np/N — 0.5, and different U. The numerical error of 
the DMRG data is of the order of the symbol size. In- 
stead of decreasing, as predicted by bosonization and in 
contrast to our findings for the other two models consid- 
ered, wo(riF] U)/wo(riF] 0) incr eases for small and mod- 
erate values of U (see Fig. 6.1). For moderate U w 2 a 
crossover to a suppression occurs at system sizes reach- 
able within DMRG. For smaller U, e.g. U — 0.5, only the 
increase can be seen. We expect that the crossover sets 
in a t mu ch larger chain length. Only the large U data of 
Fig. 6.2 display a clear suppression of the spectral weight 
for all the system sizes available and thus only for these 
data a comparison to the power-law with exponent as = 



[K 1 — l)/2 predicted by bosonization seems meaning- 
ful. K p (U, n/ ) for the Hubbard model has been calculated 
in Ref. M. For the parameters chosen we have K p (U = 
8,n/ = 0.5) = 0.62 and K p (U = 16,71/ = 0.5) = 0.56. 
Power-laws with exponents as{U = 8,71/ = 0.5) = 0.31 
and ccb{U = 16, n/ = 0.5) = 0.39 are shown as solid 



lines in Fig. 6.2. For large N the numerical data seem 
to approach these lines and we conclude that the DMRG 
results are consistent with a final power-law suppression 
of the spectral weight near the boundary which is given 
by the boundary exponent as- The asymptotic behavior 
is thus consistent with the prediction of bosonization [p"l[ 
|l2| and results obtained using Bcthe ansatz and boundary 
conformal field theory [|l3]. 

The surprising new finding is the increase of weight 
for small and moderate U and the subsequent crossover. 
This could have been expected already from the lowest 
order result Eq. ( [j.6[ ) for the spectral function. For a k- 



independent interaction U the leading logarithmic correc- 
tion to p H is given hypOf 



U 



2ttv f 



hi 



Vf^f 



(6.2) 



The prefactor of the logarithmic correction in a model 
with spin (z = 2) has thus the opposi te s ign as in the 
case of a long range interaction Eq. ( |4.2| ). As long as 
[U/(2ttvf)] In (N) -C 1 this indicates a logarithmic in- 
crease of the weight displayed in Fig. 6.1. The crossover 
scale can be estimated from this expression to be ~ 
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exp [— 2-kv f /U\. From the data it is clear that also the 
exact crossover scale 1/N C strongly depends on U. Going 
from U = 2.5 to U = 2 it roughly decreases by one or- 
der of magnitude. As a consequence of this and the fact 
that we are limited to N = 256 lattice sites the crossover 
can only be observed within a small window of interac- 
tions 2 < U < 3. Below we will show that the crossover 
scale indeed decreases exponentially in — 1/U using our 
analytical Hartree result and scaling arguments. The 1/N 
dependence of u>o(n F ; U)/wo(n F ; 0) can be translated into 
the lo dependence of Pi(uj) by multiplying by the inverse 
of the noninteracting density of states irv F . For the en- 
ergy scale at which the decrease of the weight sets in we 
thus find A = v F ir/N c . For e.g. the U = 2.5 data at quar- 
ter filling this leads to A/B « 10 -2 , where B denotes the 
bandwidth B = 4. This implies that for the given parame- 
ters the prediction of bosonization only holds for energies 
much smaller than one hundredth of the bandwidth. In 
Sec. @ we will further discuss this issue. 



1.3 



S. 1.2 

o 

g 
S 

o 

5 1.1 



■ 


a 


■ 








♦ 


♦ 


♦ 




1 ♦ 


-A 












■ 


• 


• 




X 












• 








A 




• • U=0.5 








/ 




■ U=2.0 




• 








♦ ♦ U=2.5 












A-— A U=4.0 







10" 



10" 



1/N 



Fig. 6.1. DMRG results for the spectral weight at fi and lattice 
site 1 for the Hubbard model at quarter filling. The interaction 
U is given in the legend. 



6.2 Numerical results: Hartree approximation 



Numerically the Hartree approximation (the Fock term 
vanishes) for the Hubbard model is best studied in the 
site representation. It can be done in close analogy to the 
l atti ce model of spinless fermions discussed in Eqs. ( p"^ ) to 
(5.4). For the Hubbard model the hopping is not modified 
by the interaction and the one-particle problem which re- 
mains to be solved for an electron of spin species s is the 
dynamics in an external Hartree potential generated by 
(%,-s)o- The Hartree potential shows Friedel oscillations, 
i.e. oscillates as cos (2k F j) and slowly decays as 1/j. In the 
continuum limit external potentials of this form are called 
Wigner-von Neumann type potentials pl[. From scattering 
theory it is known that the eigenvalues and eigenstates of 
Hamilton operators with oscillating and slowly decaying 
(~ 1/x) external potentials have unusual properties f3lf . 




Fig. 6.2. The same as in Fig. 



but for larger U. The solid 



lines are power-laws with exponent as = (K f 



l)/2. 



Numerically the spectral weight of the lattice model can 
be calculated along the lines discussed in the last section. 
In Figs. 2 and 3 of Ref. we have presented numerical 
data for (uj) as a function of uj, different U , and system 
sizes of up to N = 10 6 lattice sites. We have restricted our- 
selves to energies close to the chemical potential. For small 
and moderate U the spectra display qualitatively the same 
behavior as the DMRG data discussed above. For energies 
approaching the chemical potential from below the weight 
first increases following a power-law. From the numerical 
data the exponent can be determined with high accuracy 
to be —U I (2i rv F ) . This shows that the logarithmic di- 
vergence Eq. ( |6.2[ ) obtained in the leading order solution 
of the Hartree problem can be resummed to produce a 
power-law in an energy range close to the chemical poten- 
tial. The numerical Hartree data then display a crossover 
on a scale which decreases exponentially in —1/U, and 
a subsequent power-law suppression. From the numerical 
data the exponent is found to be 



H 



U/(2irv F ). 



(6.3) 



has to be compared to the leading behavior of the 
exact boundary exponent olb = {Kp 1 — l)/2< It is given 
by a B = U/(4nv F ) + O ([U/2ttv f } 2 ) |§ which is one-half 
of atg. This kind of discrepancy between exponents ob- 
tained in perturbation theory and the leading behavior of 
exact exponents is known from the Hubbard model with 
PBCj33|. It occurs because the scaling of coupling con- 
stants (of the backscattering contribution in case of PBC) 
is not taken into account in simple perturbation theory. 
Thus the coupling constant relevant for exponents is over- 
estimated in a naive perturbative treatment. For the same 
reason DMRG and Hartree data only show qualitatively 
agreement even in the limit of small coupling (see Fig. 4 
of Ref. ^jj). This has to be contrasted to the above re- 
sults for the TL model and the lattice model of spinless 
fermions. In analogy to PBC the couplings in these two 
models do n ot scale H. 

Fig. |6.3| shows pf{oj) as a function of uj for all ener- 
gies within the band, N = 2000, U = 5, and two fill- 



V. Meden et al.: Luttinger liquids with boundaries: Power-laws and energy scales 



11 



ing factors nf which are related by nj — > 2 — nf. The 
individual weights have been connected to a continuous 
line. Since the hopping amplitude is not renormalized by 
the interaction the total bandwidth is equal to the non- 
interacting one B = 4. Besides the crossover behavior for 
energies close to the chemical potential, i.e. to — 0, the 
solid line with nf = 0.8 shows a symmetric suppression 
of weight around to = —2s (hp) > 0, i.e. at energies which 
are unoccupied in the ground state. On the Hartree level 
the spectral function shows two important symmetries. In 
the next subsection the reason for both, the supression at 
to = —2e(kp) and the symmetries will become clear. The 
spectral function for fixed filling and repulsive interaction 
U > can be mapped onto the one for attractive interac- 
tion with the same absolut value — \U\ and the same filling 
by taking the mirror image at to = —e(kp). In this way 
the weight for U > around to — —2e(kp) is mapped onto 
the one for — \U\ around to = 0, i.e. the chemical potential. 
Here we do not consider the spectral function for negative 
U any further since the Hubbard model with attractive 
interaction is n ot a LL[^,|^]. The second symmetry is ob- 
vious from Fig. |6.3| . The spectral function for fixed U and 
filling rif < 1 can be mapped onto the one with filling 
2 — rif by taking the mirror image at to = 0. Note that 
the weight for U > and energies around the chemical 
potential shows a strong asymmetry. It implies that for 
rif < 1 the increase preceding the final suppression of the 
weight of energies which are occupied in the ground state 
is more pronounced compared to the analogous weight for 
fillings nf > 1. This has to be contrasted to bosonization 
which always gives symmetric behavior around to = 0. 

Next we will present an analytical discussion of the 
spectral weight within the Hartree approximation. We will 
be able to analytically determine the exponents of the 
power-law increase, the subsequent decrease and the cross- 
over scale A(U,rif) discussed above. 




Fig. 6.3. Spectral function of the Hubbard model at site j = 1, 
for TV = 2000, (7 = 5, and two different filling factors n/ in the 
Hartree approximation. 



6.3 Analytical results 

For an analytical discussion of the Hartree approximation 
it is preferable to work in fc-space. Sim ilar to th e case of 
a general continuum model Eqs. (2__3) to (2/7) we thus 
first calculate the matrix elements of the local Hubbard 
interaction in the noninteracting eigenstates 



<fn(J) 



N + l 



sin (fc„j) 



(6.4) 



with k n = mr/(N + 1), n G {1,2,..., N}. They can be 
written similar to Eq. ( J2.E| ) with 



'(q m ,q m ') = S s 



U 



N 



N 



cos (q m j) cos (qw j) 



4^ 



m,—m' } ' 



(6.5) 



and the (N + l)-periodic Kronecker 5, 5^ m , — 1 for m = 
m'+ l(N +l), I € Z and zero otherwise. A comparison with 
Eq. (2/7) shows that the correction term g(q,q') vanishes 
for a local interaction. From E q. (|6.5| ) the Hartree self- 
energy can be calculated. In Fig. 6.4 the self-energy matrix 



Ef = Ef - 5^1 



(6.6) 



with 5[i = U(n F + l/2)/(N + 1), is sketched for a filling- 
factor rif < 1. S imilar to the TL model the function f([n+ 
n']/2) Eq. (3/2) enters the expression for the self-energy 
and only matrix elements with even n + n' do contribute. 



In the hatched area and for even '/ 



the matrix elements 



are given by — U/[2(N+ 1)]. They vanish outside this area. 
For half filling nf — 1 all the off-diagonal elements of 
the self-energy matrix vanish and the spectral function is 
given by the noninteracti ng o ne. For fillings rif > 1, E^ is 
given as sketched in Fig. 6.4 but with np —> N — np and 
U — » —U. This leads to the symmetries of the spectral 
weight discussed in the last subsection. From now on we 
will restrict ourselves to fillings rif < 1. The weight for 
rif > 1 can be constructed using the discussed symmetry. 
From the self-energy the Green's function [G(to)] nn , can 
be calculated following Eq. ( |3.3| ). 



1 n F 2n 




-U/(2N+2) 



N-2n F N-n F 

Fig. 6.4. The self-energy matrix for the Hubbard model 
with n/ < 1. For details see the text. 



A lowest order inversion of [G H ] leads to Eq. (|J). 
The calculation shows that the nonanalytic behavior of the 
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spectral function close to the chemical potential, i.e. close 
to u) = 0, comes from the step through the point (wp, np) 



in the upper left part of the self-energy matrix Fig. 3.4. It 
is given by the line n' = 2np + 1 — n. The sharp step in 
the lower right half, given by n' = 2(N + 1) — 2n F — 1 — n, 
com es fr om umklapp processes where m = ml + N + 1 in 
Eq. ( J6 - 5| ) . In the spectral function it leads to the symmetric 
suppression of the weight close to the unoccupied energy 
lo = —2e(kp) already discussed in the last subsection (see 
the solid line in Fig. |3.3] ). Numerically it can be shown that 
the suppression is given by a power-law with exponent 
U/{2irv F ). 



6.3.1 Effective model 

As we are only interested in the nonanalytic behavior of 
the spectral weight for energies close to the chemical po- 
tential the Hartree problem can be simplified. This leads 
to an effective low-energy model. Due to the structure of 
the self-energy matrix the Hartree problem separates into 
two equivalent problems. One for even n and n! and the 
other one for odd n and n' . In both problems the vanish- 
ing matrix elements in the hatched area of Fig. |6.4| are 
left out, but the level spacing is doubled. Equivalently 
one can work with a "coarse grained" version of . It 
has nonvanishing matrix elements for all n and n' in the 
hatched area of Fig. x4 , but they have only half the size 
of the original ones. As only the sharp step in the upper 
left part is important for the nonanalytic behavior close 
to yU the matrix can be replaced by a matrix V with 
elements 



u 



4(iv + i; 



0(n + n' - 2n F ). (6.7) 



Next the space of one-particle states is changed. Without 
modifying the low-energy properties it can be reduced to 
n € {1, 2, . . . 2np}. Furthermore the one-particle disper- 
sion £(fc) can be linearized around kp 



£(fc — kp) — vp(k — kp). 



(6.8) 



In a last step the quantum numbers n are shifted by np. 
We then have to determine the local spectral function for 
a scattering problem with the Hamiltonian 



H 



-—np + 1 



v F k n a\,a n + 



E 

n,n' — — hf-\-1 



V n ,n' < a „(6.9) 



In order to obtain a continuous spectral function we take 
the thermodynamic limit. This is accomplished by switch- 
ing to unperturbed one-particle states 



and using 



0fe(i) = V - sin [( fc « + k F )j]. 



(6.10) 



(6.11) 



In these states the scattering potential is given by 

U 



V fe „ 



6.3.2 Scattering theory 



4tt 



G(k + k') 



(6.12) 



As usual in scattering theory |Q we define the off-shell T 
matrix T(uj) which obeys the Lippmann-Schwinger equa- 
tion 



T{uj) = V + VG {u)T{u), 

with the bare retarded Green's function 

1 



(6.13) 

S(k-k'). (6.14) 



uj — vpk + iO 

The loca l spectral function can then be calculated from 
Eq. (O) and 



G(uj) = G (u) + Go(u)T(u)G (u). 



(6.15) 



For the present problem it turns out to be advantageous to 
express the Green's function in terms of the so called K (or 
Heitler) matrix [j34| which obeys the Lippmann-Schwinger 
equation 

K(u) = V +VG$(u)K(u), 



with 



V 



k,k> 



vpk 



d(k-k'). 



(6.16) 



(6.17) 



Here V denotes the principal value. T and K are related 

by 



\K(uj)] h - \K(lu)]~ ,, 

Vp + ITT [K{u))\ a a 



with Cj = Lo/vp. Using Eqs. ( 6.15 ) and ( 3-4 ) this leads to 

p f( u )= p °(u;) -^M (6.19) 
J J 1 + ir z g(oj) 

wher e the noninteracting spectral weight is given by Eq. 

(EED 



p» = 



1 — cos (2 [kp + uj/vp] j) 



TTVp 



(6.20) 



Close to the boundary and for \uj\ -C v F k F , p®(ui) can be 
replaced by its value at lo — 0. The functions g and hj are 
given by 



g(u) = [K(w)]l a /4 



and 



hj{uj) 



-k F <Ps>U) u-vpk 



(6.21) 



(6.22) 



with (f>k(j) as in E q. ( | 6.10 ). The leading contribution to 
the integral in Eq. ( |6.22j ) comes from the pole at to — vpk. 
For small j, i.e. lattice sites close to the boundary, small k 
and small ui, <PkU)/ < P&U) depends only weakly on k and 
can thus be replaced by unity. 
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6.3.3 Spectral function 



7 Discussion and Summary 



In the Appendix we partly solve the L ip pman n- S chwinger 
equation for the off-shell K matrix Eq. (6.16) and present 
arguments which strongly suggest that the leading behav- 
ior of g(co) and /ii(w) for small \uj\ is given by power-laws 



U 



and 



Sirvp 



2uj 



vpk} 



sign(w) 



2lo 



(6.23) 



(6.24) 



The boun dary exponent a.g is defined in Eq. (6J3). Using 
Eq. ( 6.1 9| ) for \uj\ -C Vpkp the spectral function at the 
boundary site is thus given by 



Vjrkf 



l 



VF^F 



sign(w) 



r.(6.25) 



The way the effective model was constructed this analyt- 
ical result captures the leading small \u>\ behavior of the 
local spectral function of the Hubbard model within the 
Hartree approximation. All t he features found in our nu- 
merical calculations Sec. 



3.2 



and Ref. |29| are confirmed 
by Eq. (|6.25| ). For U > and \u\ -^Owc first find a power- 
law increase of the we ight with exponent -af given by 
the numerator of Eq. ( |6.25 ). For even smaller \u>\ the de- 
nominator becomes impor tant and a power-law decrease 
with ctg sets in. From Eq. (6.25) the crossover scale A can 
be calculated analytically 



Vpkp 



exp 



TTVp 1 + [U/(8v F )f 



U 



(6.26) 



Up to logarithmic corrections it is exponentially small in 
— 1/U as already expected from the numerica l dat a. The 
sign(w) function in the denominator of Eq. (|6.25 ) is re- 
sponsible for the asym metry of the spectral weight around 
oj = discussed in Sec. 3.2 

The analytical result Eq. (3.25) makes clear that it 
was essential to dete rmine the leading small M behavior 



of the K matrix Eq. ( |6.16 ), respectively of g(u>) Eq. ( 6.21 ) 



and hj{uj) Eq. ( 6.22| ), instead of the leading behavior of 



the T matrix: The latter one gives the final power-law 
decrease, but does not capture the power-law increase and 
the crossover. 

For U < only the numerator of Eq. ( |6.25 ) is relevant 
and leads to the symmetric power-law suppression of the 
weight close to /i with expone nt \ U\/ (2nvp). This confirms 
the numerical results of Sec. |6.2| . 

We note in passing that by calculating the local spec- 
tral function of the ID Hubbard model in the Hartree ap- 
proximation, we have determined the local spectral func- 
tion for the one-particle scattering problem of an electron 
in an external cos (2k f x)/x potential. The Jost function 
and scattering matrix for this problem have recently been 
discussed in the literature [fill. 



Our study of the one-particle properties of LL's with open 
boundaries covers different aspects of the low-energy physics 
of ID electrons with repulsive interaction. We have care- 
fully reinvestigated the bosonization approach to ID mod- 
els with OBC, especially focusing on the subtleties, which 
go beyond the ones known from bosonization of models 
with PBC. As no RG study exists which takes all the 
scattering processes into account we have restricted our- 
selves to interactions which are long range in real space 
and the low-energy subspace of the Fock space applying 
Tomonagas original ideag to models with OBC. We were 
then able to confirm the result for the boundary exponent 
cub which characterizes the suppression of spectral weight 
close to the boundary and the chemical potential |ll| . olb = 
(Kp 1 — 1)/ z can indeed be expressed in terms of the bulk 
LL parameter K p of the Tomonaga-Luttinger model. Led 
by this observation and analogy to PBC several authors 
have generalized the above results to all models of LL's[|il"[ 
|l2"| . Using the numerically exact DMRG algorithm we in- 
vestigated whether this is a legitimate generalization for 
two models with a short range interaction: The lattice 
model of spinlcss fermions with nearest neighbor interac- 
tion and the ID Hubbard model jD|[l4|. In both cases we 
were able to explicitly verify that the final suppression 
of the spectral weight at small energies is consistent with 
the prediction of bosonization, although for the Hubbard 
model only at energies surprisingly close to the chemical 
potential, respectively for very long chains. 

Besides that our study has revealed two interesting 
new results. Firstly we found that many aspects of the 
influence a boundary has on the low-energy one-particle 
spectra can be understood within the Hartree-Fock ap- 
proximation for the self-energy. Using numerical and an- 
alytical techniques we were able to explain how the non- 
trivial structure of the self-energy matrix leads to power- 
law behavior of the Hartree-Fock spectral function close 
to the boundary. For the two models with dominant for- 
ward scattering (TL model and lattice model of spinlcss 
fermions) quantitative agreement between the exact spec- 
tral function and the approximated one can be reached 
for small interactions. In both cases the exact boundary 
exponent olb and a§ F agree up to leading order in the 
interaction and the energy range over which power-law 
behavior can be observed is large compared to the one 
observed in the Hubbard model. 

The second surprising finding is the crossover behav- 
ior of the spectral function of the Hubbard model. In both 
the numerically exact DMRG calculation and the Hartree 
approximation the spectral weight for small and moderate 
U initially increases for energies approaching the chemical 
potential. Only below a crossover scale A the power-law 
suppression predicted by bosonization sets in. The DMRG 
data reveal that A strongly decreases with decreasing in- 
teraction U. Within the Hartree approximation we were 
able to show analytically that up to small corrections 
A — Aq exp (— Uq/U), with Uq — irvp and Aq — Vpkp. 
Due to the scaling of coupling constants not captured in 
the Hartree approximation only qualitative agreement be- 
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tween the DMRG and Hartree spectral weight can be 
reached. The scaling is responsible for the fact that the 
Hartree boundary exponent a 1 ^ is twice as large as the 
leading behavior of the exact exponent. Furthermore the 
crossover scale is underestimated by the Hartree approx- 
imation (see Fig. 4 of Ref. |2l|). Nonetheless we believe 
that the exact crossover scale shows an exponential — l/U 
dependence, with a modified Uq and Aq. Within a per- 
turbative RG calculation similar to "g-ology" we expect 
that irrelevant couplings scale to zero as 1/ |ln(yl/ylo)| if 
A — > 0. Here A denotes the momentum scale up to which 
degrees of freedom have been integrated out and A is the 
bare momentum cutoff (both measured relative to k F )- 
Thus irrelevant couplings scale down on exponentially (in 
— l/U) small scales. This implies that only Aq and Uo of 
the exponential dependence of A found in the Hartree ap- 
proximation are modified. A similar result can be obtained 
if one starts at the LL fixed point and takes into account 
the anomalous dimension of the fermion field. Close to 
the LL fixed point the leading irrelevant couplings are ex- 
pected to scale to zero as (A/A )^ u \ with j(U) ~ U. 
This again leads to an exponential scale on which irrele- 
vant couplings scale down. Our observation has dramatic 
implications for the energy range over which the prediction 
of bosonization is valid in the present context: It vanishes 
exponentially in —l/U. 



The perturbative result Eq. (|3.6| ) for the spectral weight 
indicates that the crossover behavior can be found in all 
models with dominant total backscattering zV(2k F ) > 
V(0). 

Over the past ten years several groups have attempted 
to find LL behavior in a variety of systems which be- 
have as quasi one-dimensional conductors using angular 
integrated and angular resolved photoemission spectros- 
copy jlj|20,2l], ^,^3). Here we do not want to give a com- 
plete account of all the experiments and the subtleties of 
their interpretation! 19 and focus on one aspect, which 
we believe has to be reconsidered in the light of our re- 
sults. The angular integrated spectra of all the possible 
candidates for LL behavior show a suppression of spectral 
weight close to the chemical potential p9|. If the data are 
interpreted in terms of the LL picture surprisingly large 
exponents a of the order of 1 are required. Within bulk LL 
theory this implies -Kp's smaller then the once for which 
umklapp scattering becomes relevant taking the filling fac- 
tors of the experimental systems. This leads to a contra- 
diction since the data do not show signs of a gap. Further- 
more in a recent angular resolved measurement |20|| the 
authors were unable to detect any dispersion of the peaks 
measured. Motivated by the theoretical expectation that 
a LL chain with impurities will, at low energies, behave 
as if it is cut into several disconnected pieces with open 
ends||,[l(| both the missing dispersion and the large expo- 
nent have been interpreted in the light of this expectation. 
If the electrons come from localized regions of the chain, 
they do not display any dispersion. For electrons coming 
from positions close to an open boundary olb would be 
the relevant exponent. As we have seen earlier olb > ol. 
This would resolve the apparent contradiction mentioned 



above, since a as of the order of 1 can still lead to a K p 
larger then the critical value at which umklapp scattering 
leads to a gapped system. Thus it has been speculated [fl2[ 
|l4|,[l9],|2(||23| that data have to be compared with calcula- 
tions for LL's with open ends. As we have demonstrated 
above in these systems the energy range over which the 
LL exponent olb can be observed might be very small and 
taking thermal and experimental broadening into account 
even masked by a pronounced peak at higher energies. At 
the lower end this energy range is furthermore cut off by 
non LL effects as e.g. interchain hopping ]36|. For a va- 
riety of reasons our results can certainly not directly be 
applied to experimental spectra. One is that for the Hub- 
bard model studied boundary exponents olb ~ 1 cannot 
be reached. What would thus be very desirable is a cal- 
culation of the boundary spectral weight in a microscopic 
lattice model with an interaction which is of longer range 
leading to boundary exponents of the order of one. We 
conclude that a convincing interpretation of the data in 
terms of boundary effects is still missing and requires fur- 
ther theoretical investigations. 
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ICTP in Trieste. He would like to thank YU, Lu for the hos- 
pitality in his group. Part of the calculations were carried out 
on a T3E of the Forschungszentrum Jiilich. 



Appendix 

In this appendix we will determine the leading small \lo\ 
behavior of the fu nctions g{to) and h\{uj) defined in Eqs. 
( 3.21 ) and ( 6.22 ) fr om t he solution of the Lippmann- 
Schwinger equation (3.16). For the dimensionless K ma- 
trix K(uj) = K{uj)/vf it reads 



k,k' 



Vk,k> + 



dp- 



p,k' 



UJ — p 



with 



V k , k >=V kik ,/v F = -UG(k + k') ) 



, (Al) 



(A2) 



wher e U — U/(4otvf), and Co = to/vp. By differentiating 
Eq. ( Al) with respect to A; a system of coupled differential 
equations for 



T x {k) = T x {k\k' ,w) = K(io) 



k,k' 



V fcl . 



and 



T 2 {k) =F 2 (k;k',u) = T 1 (-k;k',Lo) 
can be derived 

m) = u> e{k '- k) -v T2{k) 



T' 2 {k) 



£r2 0(k' + k ) , -^i(A) 



G) + k ' 
U- 



(A3) 

(A4) 

(A5) 
(A6) 
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The boundary conditions are Ti(—k F ) = J-2(kp) = 0. 

We did not succeed in analytically solving this system 
of differential equations for arbitrary k, k', and small \u>\. 
To determine the leading small \u>\ behavior of g{uj) and 
h\(uj) we fortunately only need T\(k; k' , uj) for cer tain spe- 
cial c ombinations of the arguments. From Eqs. ( |A3[ ) and 
( |6.21 ) it follows that g(uj) is given by 

1 2 



g(uj) = Ti(uj;uj,uj) + Vz,z 



(A7) 



The differe ntial equation (A6) can be used to express 
hx(ui) Eq. (|6.22D in terms of ^(fej?; w,w)[§5| 



dk 



Fi(k;uj,uj) + V k ,u> 



-1 2 



UJ — k 
-i 2 



1 



-yFi(k F \ti,u) 



(A8) 



We thus only have to evaluate JFi (fcp ; tD, uj) and (w; w, uj) 
in the limit |w| <C vpkp. 

In a first step Eqs. ( |A5| ) and ( |A6| ) are solved for uj = 0. 
In this case a new variable t{k) = In |fc| can be introduced 
and the equations can be decoupled by introducing T± = 
T\±Ti. The rest of the calculation is straightforward and 
leads for \k\ > \k'\ to 



F 1 (k;k',uj = 0) 

U_ 

^ ~2 



U0(k 
u 



+ k') 
sign(/c) 



k 


u~ 




k' 


kp 






k F 



(A9) 



For \k'\ > \k\, k and k' have to be interchanged. I f we 
set k = kp, k' = uj, and assume \uj\ <C vpkp, Eq. (A9) 
simplifies to 



F-L(k F ;u),0) = U-U 



U 2 In 



vpkj 



-V 



vpkp 



U\ 2 



vpkj 



(A10) 



The first few terms in the expansion in powers of U can 
easily be obta ined by iterating the Lippmann-Schwinger 
equation ( |6 .16 ) for the K matrix. Unfortunately the quan- 
tity presented in Eq. ( A10 ) is not quite what is needed to 
calculate hi(ui). An exact solution is difficult, but it is 
straightforward to calculate the first terms in the U ex- 



pansion for K(uj) 



F\(kp; Co, to) = U In 



For <§; v F kp this yields 



vpkf 



U\ 2 



2uj 



vpkf 



If we now assu me that this expansion can be resummed 
as in Eq. ( AlCj ) we obtain 



/ii(cj) 



2u 



v F kp 



-2U 



(All) 



For k = uj and k' = uj Eq. (A9) gives 



TiiCj-Cj,®) = U0(oj) 



= U z In 



U 
1 



UJ 



v F k} 



-2U 



sign(w) 



vpkf 



U 3 In 2 



Iteration of Eq. ( 6.16 ) for K(uj 
2oj 



Ti(uj; uj,uj) — U In 



v F k F 



-t/ 3 ln 2 



v F k F 
leads to 
2uj 



(A12) 



v F k F 



If we again assume tha t the logarithm can be summed to 
a power-law as in Eq. ( A12) we obtain the result for g(uj) 
presented in Eq. (|6.23 ). 
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